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Abstract 

The light top-squark ti may be the lightest squark and its hfetime may be 'long enough' in a kind 
of SUSY models which have not been ruled out yet experimentally, so colorless 'supersymmetric 
' hadrons (superhadrons)' (tiq) {q is a quark except t-quark) may be formed as long as the light 

\D 

' top-squark ti can be produced. Fragmentation function of ti to heavy 'supersymmetric hadrons 



X 



(superhadrons)' {tiQ) (Q = c or h) and the hadronic production of the superhadrons are investigated 
Q-i! quantitatively. The fragmentation function is calculated precisely. Due to the difference in spin 

of the SUSY component, the asymptotic behavior of the fragmentation function is different from 
those of the existent ones. The fragmentation function is also applied to compute the production 
of heavy superhadrons at hadronic colliders Tevatron and LHC under the so-called fragmentation 
approach. The resultant cross-section for the heavy superhadrons is too small to observe at Tevatron, 
but great enough at LHC, even when all the relevant parameters in the SUSY models are taken 
within the favored region for the heavy superhadrons. The production of 'light superhadrons' {tiq) 
{q = u, d, s) is also roughly estimated. It is pointed out that the production cross-sections of the 
light superhadrons {tiq) may be much greater than those of the heavy superhadrons, so that even at 
Tevatron the light superhadrons may be produced in great quantities. 
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I. INTRODUCTION 



Supersymmetry (SUSY) is one of the most appealing extensions of the Standard Model 

(SM)ll|, mm 

. Whereas without knowing the realistic SUSY breaking mechanism, even in 
the minimum supersymmetry extension of the Standard Model (MSSM), there are too many 
parameters in the SUSY sector, which need to be fixed via experimental measurements. If 
the MSSM is rooted in the 'minimum supergravity model' (mSUGRA), the numbers of the 
independent parameters can be deducted a lot, but there are still many un-fixed parameters s], 
y,0]. Therefore, the spectrum of the SUSY sector in SUSY models still is a quite open problem^. 

For some kinds of SUSY models and by choosing un-fixed parameters from the region which 
has not been ruled out yet, it is not difficult to realize a general feature of the two mass 
eigenstates ti and ^2 for the SUSY partners of top-quark (top-squark and tRj such as that 
the comparatively lighter one ti is the lightest squark, and the Fi's lifetime is so 'long', that 
its width is less than Aqcd^, a, 11, 12, 13, 14, 15]. In this case, ti may form various 
colorless hadrons i.e. the superhadrons by QCD interaction, which consist of ti and q (here 
q = u,d, c, s, b). On the other hand, the direct experimental searching for the SUSY partners 
may only set a lower bound of ti mass as mj^ > 100 Gev[l6|,0,Q2 

(even lower as 75 GeV|15[|). 



Therefore in the paper, we would like to focus our attention on the consequences for the possible 
features of ti. Namely, we shall assume that ti, the SUSY partners of top-quark, is not very 
heavy e.g. mj^ ~ 120 ~ 150 GeV and has a 'quite long' lifetime Tj^ < Aqcd, that ti (after 
producing and before decaying) has chances to form colorless superhadrons (tiq)^. Moreover 
we think that the squark (anti-squark) in the superhadrons is a scalar, which is different from a 
quark in 'common' hadrons, hence with such a scalar component the study of the superhadrons 
is also very interesting from the bound state point of view. 

There was a remarkable progress in nineties of the last century in perturbative QCD (pQCD) 
for double heavy meson studies, i.e., it was realized that the fragmentation function and the 
production of a double heavy meson such as and 7]^, J/ip can be reliably computed in terms 



of pQCD and the wavefunction derived from the potential model 22 



23 



24 



25 



m 



27 



28|^ 



^ In fact, all of the available indications on the masses of the SUSY partners are abstracted from experimental 
measurements and/or astro-observations under assumptions (not direct measurements), so one should consider 
them only as references. 

^ For a summary see |16| . 

^ This kind of superhadrons are bound states of a quark (anti-quark) and an anti-squark (squark), or two 
gluinos, or two quarks (antiquarks) and a squark (anti-squark), etc. All of them are colorless and bound via 
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2li 



strong interaction (QCD) 

^ Exactly to say, here we mean color-singlet mechanism only i.e. only the color-singlet component of the con- 



and further progress in formulating the problem under the framework of the effective theory: 



non-relativistic QCD (NRQC D') |29l.l30[| was made a coup 



before Re: 



meson 
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28| there were papers 



e years late. We should note here that 



31 



32| . The inclusive production of a 



32| and the fragmentation functions of a parton into a meson 32| were calculated. 



But the authors of 



32| | precisely claimed that their calculations might be extended to the cases 



for the heavy-light mesons, such as D^B etc.. In fact, the claim is incorrect and misleading in 
the key point on the theoretical calculability of the production and the fragmentation functions^. 
Since heavy-light mesons contain a light quark, and the light quark creation involving in the 
fragmentation function is non-perturbative, so it cannot further be factorized out a hard factor 
as that in the case of the double heavy mesons. It is just the reason why, similar to the case for 
double heavy mesons, we expect that only the 'heavy superhadrons' {tiQ) (but not the 'light 
superhadrons' (tiq)), where Q (g) denotes a heavy quark, c or 6 (a light qluark, u or d or s), 
and their inclusive production and fragmentation functions may be calculated reliably. The 
fragmentation functions to the 'heavy superhadrons' {tiQ) can be simply attributed to a wave 
function of potential model and a pQCD calculable factor as in the case of the double heavy 
mesons. 

For the 'light superhadrons' (tiq) where q indicates a light quark: m or or s, due to 
the non-perturbative nature for producing the involved light quark q, the 'story' about the 
calculation of the fragmentation function is very different, i.e., the fragmentation function 
cannot be attributed to a wave function and a hard factor of pQCD. Due to the non-perturbative 
QCD effects in the fragmentation function of the 'light-heavy mesons' such as B, D and etc, 
so far practically the way to obtained the fragmentation functions of the 'light-heavy meson' is 
that they are formulated in terms of theoretical considerations and parametrization first, and 
then the parameters in the formulation are fixed via experimental measurement (s). With the 
fragmentation functions the production cross-section of a 'light-heavy meson', in experiences, 
generally is greater than that of the respective double heavy meson (the q quark in 'light-heavy 
meson' is replaced by a charm quark c) by a factor 10^~^. Since there is no experimental 
observation on superhadrons, so we cannot play the same way for the 'light superhadrons' as 



cerned double heavy meson, which is the biggest in Fock space expansion, is taken into account in calculating 
the fragmentation function and production as well. While since the color-octet matrix element appearing in 
the formulation for color-octet mechanism production (fragmentation function) cannot be calculated theoret- 
ically so far, so the color-octet mechanism is not the case. 
^ To present the calculations of the fragmentation functions, we would like also to remind here that there are 
some substantial contributions from the phase-space integration which might be missed if enough care were 
not paid. In fact, as pointed out in 23 1, they were missed in ^S^l indeed. 



'light-heavy mesons' at all. Alternatively, as a magnitude order estimate, we expect that the 
fragmentation function and the production of the light superhadrons [tiq) are also greater than 
that of the respective heavy superhadrons (tiQ) by a factor 10^""^, no matter how heavy ti is, 
that is very similar to the case of a double heavy meson vs a heavy-light meson. Thus, based 
on the quantitative computation of the fragmentation function and the production of the heavy 
superhadrons, we simply extend the results of the the production to the light superhadrons at 
the end of the paper by referring the cases of the 'double heavy mesons' vs the 'light-heavy 
mesons' experientially. For convenience, later on we will denote (tiQ) as H throughout the 
paper. 

Since the spectrum and the wave function accordingly of a double heavy-quark binding 
system i.e. a system of a heavy quark and a heavy anti-quark system, (Q'Q), can be obtained 
theoretically in terms of non-relativistic potential model inspired on QCD quite well, so the 
'heavy superhadrons' H, as the double heavy system (Q'Q), may also be depicted by the non- 



re. 



2Q|, 



ativistic potential model as long as the difference in spin is taken into account carefully 
2l|. Therefore, there is no problem to obtain the wave functions of 'heavy superhadrons' H 
that appear in the fragmentation functions. 

In the literature, there are two approaches for estimating the direct production of a double 
heavy meson in NRQCD framework: the 'fragmentation approach' vs the complete 'lowest- 
order-calculation' approach. It is known that, of the two approaches, the former is much 
simpler than the latter one in computation, but the former is 'good' only in the region where the 



transverse momentum of the produced double heavy meson is large {px ^ 15GeV)[33 



34 



35|. 



The situation for the production of the heavy superhadrons is similar to the cases of the double 
heavy mesons, so of the two approaches, we adopt the fragmentation approach when estimating 
the production of the superhadrons for rough estimation. 

This paper is organized as follows: In Section II, we show how to derive the fragmentation 
function of the lightest top-squark ti to the 'heavy superhadrons' H, and try to present the 
obtained results i.e. its general features properly. In Section III, we compute the cross sections 
for hadronic production of the superhadrons at colliders Tevatron and LHC in terms of the so- 
called fragmentation approach. The Section IV is devoted to the discussions and conclusions. 



II. FRAGMENTATION FUNCTION OF THE LIGHT TOP-SQUARK h TO THE 
HEAVY SUPERHADRONS H 



In this paper we adopt the fragmentation approach to estimate the H production, and in 
the present section, we are computing the fragmentation function first, that is one of the key 
factors of the fragmentation approach. 

According to pQCD, with leading logarithmic (LL) terms being summed up, a fragmentation 



function of a 'parton' i to a heavy superhadron H is depicted by DGLAP Equation as below 
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*H 

dr 



E 



2n 



y ^ 



(1) 



where r = log{Q^ / Aq^jj), Pi_,j{x) is the splitting function. For example, the splitting function 
for the supersymmetric top-squark jsol, 40 1 reads 



p~ ~ {x) = - 



:i-x) 



x] 



6{1 



x] 



Since Eq.(IT]) is an integro-differential equation, so to have definite solution, a 'boundary 
(initial) condition' for the equation i.e. D^_^j^{z, Qo), the fragmentation function at the energy 
scale Qo ~ m^^, is needed. Now the task is to obtain the boundary condition. Fortunately, 
the boundary condition for heavy superhadron H can be derived in terms of pQCD and the 



28|. Hereafter, 



relevant wave function precisely as double heavy mesons|22|, |23|, |2j, |25|, |26|, |27, 
to simplify our notation, we shall always use D_._^^{z) instead of Dj^'^{z,Qo). 

Since the fragmentation functions are universal by definition, i.e. they are independent of 
the concrete process, so for the 'boundary (initial) condition' of the fragmentation function 
of the light top-squark ti, we would like to choose a relevant simple process to calculate the 



'boundary condition' D 



ti~-*H\ 



In order to simplify the derivation as much as possible, we 



furthermore assume a fictitious "Z" , which, except the mass, has the same properties as that 
of the physical Z boson. The fictitious "Z" has such a great mass that it may decay to a pair 
of ti and ti. 

According to pQCD factorization theorem, the differential width for the fictitious "Z" de- 
caying to H may be factorized as 

dr{''Z" + X) = r dzdf{''Z" ^h + ti,fif)D~_^j^{z,fif), (2) 



where z 



2E 



and /i/ is the energy scale for factorization. By definition, Dj^_^j^{z, ^f) is the 



fragmentation function, which represents the probability of ti fragmenting into the superhadron 
with enerffv fraction z. 
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FIG. 1: The Feynman diagrams for the fictitious particle ^^Z"(k) decaying into a superhadron H(p), 
ti(gi) and a bottom quark 6(^2) (or a charm quark c((72))- (A), (B), (C) are the Feynman diagrams 
with the hard virtual-gluon attached to the Ught top-squark ti and its anti-particle ti in different 
ways. 

To calculate D~^^g{z), let us calculate the process "Z" ^ H + X precisely. Of the lowest- 
order pQCD calculation, the Feynman diagrams for the inclusive decay of the fictitious particle 
"Z" into a (|)^ superhadron H (with mass M), "Z" H + X , are described by the three 
diagrams (A), (B) and (C) as shown in FIGO The intermediate gluon in each of the Feynman 
diagrams should ensure the production of a heavy quark- ant iquark pair, so its momentum 
squared should be bigger than {p2 + ^2)^ > ^mg ^ ^%cdi ^^us pQCD calculation and the 
factorization theorem are reliable. The corresponding amplitudes are 

with 

[k — qiY — m~ [k — piY — m~ 



and 



ju X ^ . X ^ giM(g2+P2)^ + gi/.(g2 + P2)/. </>(0) 

L = M g2 7 P , G^i^ = 9t.u J ■ T 9b = 77— 

gi- g2 + P2 2^ 



where k, pi, p2, qi and q2 are the four momenta of "Z", top-squark ti, antiquark Q (b or 
c), top anti-squark ti and quark Q {b or c) respectively, p is the four momentum of H and 
q is the relative momentum between the two constituents inside H, so we have p = Pi + P2, 
q = a2Pi—CiiP2 with ai = ^^^^ , 0^2 = m^m2 ' 'where mi and m2 are masses ofti,Q respectively. 
0(0) is the wave function at origin for the superhadron H. If neglecting the g-dependence of the 
integrand of Eq.([3]), which can be considered as the lowest term in the expansion of q for the 
integrand, all the non-nerturbative effects can be attributed into the wave function at origin 



TABLE I: The wavefunction at origin (j){0) for or (tic). Here the needed parameters are set: 
rub = 5.18 GeV and = 1.84 GeV. 



mr 


120 GeV 


150 GeV 




"(Gey)3/2" 


2.502 


2.530 






0.693 


0.695 



after doing the integration over q. e is the polarization vector for the fictitious particle "Z' 
cii = Izlcos^O-— CqSin^Ow [41] ^. For convenience, the variables: S", 



2p-k 



, y 



M = m^^+ ruQ and d 



M 



2qi-k 

Q ,., 1 y fit' 

■'e// Jeff 

are introduced. Keeping the leading term for ci^, the 

nax — J- x{l-x) ■ 



maximum and minimum values of y are y^aa; - . , y^^^ - . ^ , ^^^^^^ 

calculation the axial gauge ra^ = gi^ is adopted. Under the axial gauge, it can be found that 



Vr. 



l_a;+f!!(lz£±^. In the 



only the two amplitudes Ma and Mb-, which correspond to the first two Feynman diagrams in 
FIGiH have contributions to the fragmentation function. 

According to the factorization Eq.([2]), the fragmentation function versus z at energy scale 
Qq can be derived by dividing the differential decay width with Fq: 



Fo dz ' 



(4) 



where Fq is the decay width for the fictitious particle "Z" decaying into top-squark ti and top 
anti-squark ti. Thus the result for the fragmentation function may be presented as follows: 



(5) 



where 



4W 



1 (l-^)V 



16a,(4m2)|0(O)p 



mr: 



2al{z - A)z + al{?,aiz - 2z + 2)z + Sa^ - 6ai + 6 



(6) 



Qil-aizy 

At preset, there is no experiment data for the superhadron at all, so we adopt potential 
model with the Cornell potential (— 7 + ^) to estimate the wave function at origin, 0(0). For 
definiteness, we assume that the potential of the heavy scalar-antiquark binding system is the 
same as that of double heavy quark-antiquark systems. The relevant parameters are taken as: 



^ cii is a factor in the effective coupling "Z" —ti—ti, which will not appear in the final result of the fragmentation 
function because of cancelation between the numerator and denominator in the computation formula Eq. ^ . 



K = 0.52, a = 2.34 GeV"^ [4^, mj^ = 120 or 150 GeV, = 5.18 GeV and = 1.84 GeV. 
The corresponding wave functions at origin 0(0) for the systems (tib) and (tic) are hsted in 
TABUl 

The fragmentation function D~^_^^{z) obtained, see Eq.([5]), is just a boundary condition 
for the DALAP evolution equation Eq.([T]). Solving the DGLAP equation, one may obtain the 
fragmentation function with the energy-scale evolution to Q^. The relevant Feynman diagrams 
for the boundary condition D.^^{z) with (j = q,q,g) are of higher order in as than the 
Feynman diagrams for see FIGjTf, therefore only the case with i,j = ti shall be 

taken into account in solving Eq.([T]), so as to meet the LL approximation criterion 

We solve Eq.([T]) with the method developed by Field |43|. When ^ , according to 
the Field's method, we have the following solution: 



with 



and 



D 



D 



y 

:z,Q')=D^^j^i6,z,Q' 
dy- 



K i'^D j^{Q,z/y,Q')P^,{y) 

z y 



Pa{x) 



+W ^D~^si^/y,Q')P^^j^iy) + 0{K') 
J z y 



PaJx) = 6 



X 



1 — X 



11 



1 — X loe X X 12 18 



(7) 



D 



dy 



g{a,z,Q ) = / —Dj j^{z/y,Qo) 
Jz y ' 



log(l/)) 



(an— I) 



^ y 



D^^s{b,z,Q') 



r(a/t) 
log(y))(^''-i) 



6 



daD 



T{bK) 

'a,z,Q^) 



3 3 

\ f'dbD^^g{b,z/y,Q'^ 



6- 



3 3 



^ For D^^~{z), the relevant part of the Feynman diagrams must have one more strong couphng vertex {g — > 
titi) in as than the relevant part ti ^ H in FIGdJ and for the relevant Feynman diagrams must 

have two more strong coupling vertex q qg and g — > titi in as than those ti — )■ H. 




FIG. 2: Behavior of the fragmentation function (amphfied by a scale factor 10^ for convenience) for 
the hght top-squark ti to H (here H precisely is the S-wave (tib) superhadron). The dotted and 
dashed lines are those stand for Eq.([5]), the 'initial' fragmentation function, with m~ = 150 GeV 
and m~ =120 GeV respectively. The solid and dash-dot lines are those stand for the fragmentation 
function evolving to the energy scale Q = 2 TeV (a typical energy scale) with m^^ = 150 GeV and 
nij^ = 120 GeV respectively. 



where k = ^^^2nf ^'^si^siQl) / c^siQ"^)) , 1e is Euler constant. Furthermore, at LL level we have 
the boundary condition D^^-^{z/y,Ql) = 0, thus the solution Eq.(I7j) becomes 

D^^S{Z,Q') = ^ f'^\_^s{z/y,Q')P^j^{y) + 0{K'). (9) 
Jz y 

Numerically, it can be found that the first term for Dj^^^{z, Q"^) is much greater than the other 
terms in the right hand side of the first equation of Eq.(l9]), and then it is quite accurate to 
consider the first term only. Moreover, due to the fact that the splitting function -P^^^~ must 
be suppressed greatly O^rriQ/m- ) as ti is heavy {mj^ > 120 GeV), we may safely conclude that 
P'g^H^^^ ~ ^ when is not very great. 

To precisely see the general behavior of the fragmentation function obtained, we draw its 
curves in FIGJ21 It can be found that when the mass of the top-squark becomes heavier, the 
peak of the curve for the fragmentation function increase higher accordingly. 

Furthermore, to see the characters of the obtained fragmentation function, let us compare 
it with those for quarks {Q). The fragmentation function for an anti-quark Q into a double 



heavy meson {QQ'), e.g., a bottom anti-quark b into Be, which can be found in Refs. 22, [23, 




FIG. 3: The comparison of the two types of fragmentation functions. The upper curve (the sohd 
one) is D'^^{z), and the lower curve (the dashed one) is D'g{z). The relevant parameters are taken as 
mj^ = mi, = 5.18 GeV, rric = 1.84 GeV. The definitions and the artificial assumption are taken as 
those in text. 



2J,l25|,l26|,l27|,l28j: 



Di{z) = F-,-Mz), 



where Ft, 



27Mm2 



and 



[12\2Z - 3(Ai - A2)(l - Ai^)(2 - z)\ (1 - \^z)z 



(1-Ai^)6 ^ 

+6(1 + \2zf{\ - Xizf - 8AiA2z'(l - z) 



(10) 



where ipoiG) is wavefunction at origin of -^i = -^2 = and M = + rric- 

To highlight the difference between the two types of fragmentation functions, we remove the 
irrelevant factors Fi and Fj^ away and introduce the functions D'^^{z) and D'g{z): 

D',{z) = h{z) and D'^z) = h^{z) . 

One may see the differences clearly in asymptotic behaviors of the two kinds of fragmentation 
functions that for D'^{z) and D'g{z) they are z and 2;^ as 2; — >■ respectively; and the same (1 — 2;)^ 
as z — > 1. FIGJ3] depicts the two kinds of fragmentation functions quantitatively. In the figure, 
the function D'^^{z) is taken precisely as the fragmentation function for the h quark fragmenting 
into S-wave pseudoscalar state of a double heavy meson Be or 5*, while the function D~ {z) 
is the fragmentation function for the liffht toD-snuark U fraementinET into S-wtwe sunerhadron 



H = (tic). To contrast with the difference between these two kinds of fragmentation functions, 
in FIGj3]we have assumed m-jr^ = rrih ~ 5.18 GeV artificially. 

III. PRODUCTION OF THE SUPERHADRON AT HADRONIC COLLIDERS 

Here we are adopting the fragmentation approach to estimate the production of the su- 
perhadrons at hadronic colliders. According to the NRQCD factorization theorem, the cross 
section of if-production by collision of hadrons Hi and H2, d(yHiH2^HX^ '"^'^ factorized into 
three factors as below: 

(^^HiH2^HX = Y.J J dx2 J dzfi/H,ixi,fIf)fj/H2ix2,fif) 
ijk 

■dcTij-^kxixi, X2, z; Hf, hr) ■ D^^g{z,Hf), (12) 

where i,j and k are parton species; /i/ corresponds to the energy scale where the factorization is 
made; fiR is the renormalization energy scale for the hard subprocess; dcrij^kx{xi,X2, z; fif, ^r) 
is the cross section for the 'hard subprocess' ij — > kX; D^^^{z^ iif) is the fragmentation func- 
tion of 'parton' k to H] /j/^^-j (xi, /i/) and fj/H2{^2,fJ'f) are the parton distribution functions 
(PDFs) in collision hadrons Hi and H2 respectively. In this paper, as in most pQCD calcula- 



tions, we choose jj,f = hr +Pt 5 ^^d, as a consequence of the choice, we may further 



set D^^j^{z,Q'^) = quite safely as argued above. Therefore, k in Eq. f|T2|) 'runs over' ti only, 
i.e. k = ti. 

By naive considerations, of all the possible hard subprocesses for the production {ij 
kX , k = ti), the gluon-gluon fusion g + g ^ ti + ti and the quark-antiquark annihilation 
q + q ^ ti + ti (here q and q are light quarks) are in the same order of strong coupling as, 
so they may be the most important ones for the production at the colliders: Tevatron and 
LHC The gluon component of PDFs in small x region is the greatest, so the gluon-gluon fusion 
should be the most important one at LHC; whereas, at Tevatron, due to a comparatively 
low CM energy, the energy- momentum fraction x of the gluon parton must be big enough to 
produce the titi pair, so as the case of top-quark pair production at Tevatron probably the 
components of valance quarks, instead of the gluon, play more important role (a review of this 



point can be found in Ref. 4J]). Therefore, first of all we highlight these two subprocesses for 
the production. Note that according to Ref.j^sl the production of the top-squark pair tit2 or 
^2^1 at the hadronic colliders is small, so we do not take them into account. 



® One will see later on that the cross-section of the production decrease with pT rapidly. Thus at Tevatron and 
LHC 11 f = iia . I m,~ + n'L is not hieh enoueh that m~ can be considered as zero. 



▲ 

(^) 






FIG. 4: The lowest order Feynman diagrams for the gluon-gluon fusion subprocess g + g ^ ti + ti. 



Now let us calculate the gluon-gluon fusion subprocesses first. To the lowest order (tree 
level), there are four Feynman diagrams as shown in FIGJH The corresponding amplitudes 
read 



Mf 



M, 



c 



2 ' 



ipi - k2y - rrii- 



2 ' 



^2^j.ab _ j.ba, ik2 ' k^Y g^^ - 2k^g^''^ + 2k-y^ 

/2l'-pab I rpba 



{ki + k^f 



Md = gt{T'^' + T'-)t^{k,)e,{k2)g^r 



(13) 
(14) 

{Pi-P2)i,'€-i,{ki)ey{k2), (15) 

(16) 



where e is the polarization vector of gluon. When taking the axial gauge with a fixed four- vector 
ra, the summation of polarization vector reads 

, ^.(^^ A)e.(fc, A) = - + ^.^ • 



The differential cross section is 



dcr{gg hh] 



2^,2 



IGvrs^ 



1-2A-- 
9 



2 

mA 
1-2-^(1 

As ^ 



2 

ti_ 

As 



dt. 



(17) 



where A = (i — m~ )(u — m~ )/s^. s, i and u are Mandelstam variables of the subprocess, 



s 
i 



{ki + k2)\ 
{Pi - ki)', 



u = {pi- k2f 



(18) 
(19) 
(20) 



which satisify s + u + t = 2ra~ . 

For the quark-antiquark annihilation subprocess, to the lowest order there is only one Feyn- 
man diagram and its Feynman amplitude reads: 



M 



2rT.aa {Pl - P2Mk2hM^l) 



9iT 



(ki + k2? 



(21) 



TABLE II: Hadronic cross sections (in unit: fb) for the superhadrons (tic) and (tib) with = (^) 
The parameters appearing in the estimate are taken as those stated in text. 





LHC {VS=U. TeV) 


TEVATRON (\/S = 1.96 TeV) 


Constituents 


subprocess gg 


subprocess qq 


subprocess gg 


subprocess qq 


mr = 120 

ti 

GeV 


(tic) 


114.51 


0.36469 


0.26975 


1.2E-3 


m 


30.489 


0.10374 


0.0696 


3.E-4 


mr = 150 
GeV 


(tic) 


42.176 


0.14591 


0.0537 


2.E-4 


m 


11.812 


0.0431 


0.0142 


7.E-5 



where ki and k2 are the four momenta for the quark and antiquark respectively. The differential 
cross section is obtained as below 

iral [s^ - 4sm|^ - (i - u)"^] 
daiqq^hh) = dt, (22) 

where the Mandelstam variables are defined as Eqs. ( |T8|19ll20p . To calculate the production via 
the subprocess of quark-antiquark annihilation, as stated above, we are interested in seeing the 
valance quarks' contributions, especially, the production at Tevatron, so here we are considering 
the contributions only from the light quarks in PDFs to the production. In nucleons only the 
light quarks u and d may be their valance quarks. 

The total hadronic cross section is calculated via the two subprocesses in terms of the 
factorization formulation Eq. fll2l) and with the help of the fragmentation function Eq.(l5]). The 
differential cross-section for the gluon-gluon fusion is in terms of Eq. lITTI) and the differential 
cross-section for the quark-antiquark annihilation is in terms of Eq. fl22|) . Since the present 
calculations are at the lowest order only, so the version CTEQ6L j46| for the parton distribution 
functions (PDFs) is taken, and for definiteness, we take rrib = 5.18 GeV, rric = 1.84 GeV 
and assume two possible values for mj^: mj^ = 120 or 150 GeV. In addition, the parameter 
Aqcd in the running coupling constant ag is taken as 0.216 GeV. The energy scale for the 
QCD factorization formulas is chosen as the 'transverse mass' of the produced superhadron: 

The total hadronic cross sections obtained at Tevatron and LHC are in TABLE [III From 
the table one may see that the cross section for hadronic production of the superhadron H at 
Tevatron is much smaller than that at LHC (almost by three order of magnitude), so when ti 
possesses the behaviors as assumed here and considering the final possible integrated luminosity, 
it is hopeless to observe H at Tevatron, but it is hopeful at LHC. TABLE [III also shows that 
the hadronic cross sections of the sunerhadron H at Tevatron and at LHC decrease as the mass 




FIG. 5: The distributions of the transverse momentum Pt (left figure) and rapidity y (right figure) 
for the superhadron H produced via gluon-gluon fusion at LHC. For Pt distribution, a rapidity 
cut \y\ < 1.5 is made. The upper one of the two dash hues corresponds to the distribution for 
the superhadron (tic) production with mj^ = 120 GeV being assumed, the lower one to that for 
the superhadron (tib) production with m~ = 120 GeV; The upper one of the two solid lines to 
the distribution for the superhadron (tic) production with mj^ = 150 GeV, the lower one to the 
distribution for the superhadron (tib) production with rrir = 150 GeV. 




p^(GeV) y 

FIG. 6: The distributions of the transverse momentum Pt (left figure) and rapidity y (right figure) 
for the superhadron H produced via quark-antiquark annihilation at LHC. For Pt distribution, a 
rapidity cut \y\ < 1.5 is made. The upper one of the two dash lines corresponds to the distribution 
for the superhadron (tic) production with m~ = 120 GeV being assumed, the lower one to that 
for the superhadron (tib) production with mj^ = 120 GeV; The upper one of the two solid lines to 
the distribution for the superhadron (tic) production with mr = 150 GeV, the lower one to the 
distribution for the superhadron (tib) production with rrij = 150 GeV. 
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FIG. 7: The Pt distributions of the superhadrons H produced at Tevatron via the gluon-gluon fusion 
(left figure) and the quark-antiquark annihilation (right figure) respectively. For the distributions, a 
rapidity cut \y\ < 0.6 is made. The upper one of the two dash lines corresponds to the distribution 
for the superhadron (tic) production with mr = 120 GeV, the lower one to the distribution for the 
superhadron (tib) production with mj^ = 120 GeV; The upper one of the solid lines to the distribution 
for the superhadron (tic) production with mj^ = 150 GeV, the lower one to the distribution for the 
superhadron (tib) production with mj^ = 150 GeV. 

of the light scalar top quark mj^ increasing. Moreover the cross sections for the superhadron 
production via gluon-gluon fusion are much larger than those via the annihilation at Tevatron 
and at LHC both. Hence, the contribution from quark-antiquark annihilation can be ignored 
in comparison to the dominant contribution from the gluon-gluon fusion. 

To present more features of the production, we also draw curves to show the distributions 
of the produced superhadron. The differential cross-sections vs the transverse momentum pt 
and the rapidity y of the produced superhadron via gluon-gluon fusion at LHC are drawn in 
FIG|5l while those via the quark-antiquark annihilation are drawn in FIGJ61 The distributions 
of transverse momentum px for the superhadron production at Tevatron via gluon-gluon fusion 
and quark-antiquark annihilation are drawn in FIGJ7] respectively. 

IV. DISCUSSION AND CONCLUSIONS 

Here the fragmentation function of the light top-squark ti to heavy superhadrons (tic) and 
(tib) is reliably computed as that in the case of a heavy quark to a double heavy meson. To see 
the characteristics of the fragmentation function, comparisons of the obtained fragmentation 
function for the light top-squark with those for heavy quarks are made by drawing curves with 



suitable parameters in FIG|3l When z approaches to zero, the fragmentation functions for 
the top-squark (a scalar particle) approach to zero as 2;^, instead of those for a heavy quark 
(a fermion particle) which behave as z] and both of them have a similar asymptotic behavior 
when z approaches to 1. 

Using the obtained fragmentation function of the superhadron H (either t\c or ti6) and 
under the fragmentation approach up-to leading logarithm, the cross sections and Pt (y) dis- 
tributions for H have been computed at the energies of Tevatron and LHC In the computation, 
the gluon-gluon fusion and light quark- ant iquark annihilation as the hard subprocess for the 
hadronic production of the superhadron H are taken into account precisely. When calculating 
Pt distributions, different rapidity cuts are taken, i.e. \y\ < 1.5 at LHC and \y\ < 0.6 at 
Tevatron. From the cross sections and Pt (y) distributions, one may conclude that one cannot 
collect enough events for observing the superhadron at the hadronic collider Tevatron, even if 
the parameters of the supersymmetric model are in a very favored region. In contrary, enough 
events for experimental observation of the superhadrons can be produced (collected) without 
difficulty at the forthcoming collider LHC. Namely if the coming 'new physics' is supersym- 
metric and the parameters are in the favored region of the concerned superhadrons and allowed 
by all kinds of the existent experimental observations, Tevatron is not a good 'laboratory' to 
observe the possible superhadron(s), while LHC may be a good one. Moreover, it can be found 
from TABLE [Tl] that the production cross-section via g + g ^ ti + ti is much smaller than 
that via gluon-gluon fusion subprocess at Tevatron and LHC. To compare with the top-quark 
production, let us note here that the quark-antiquark annihilation for the top-quark produc- 
tion is comparable to the gluon-gluon fusion at LHC 47|, |48|, |49|, l50|, l5l|], while for the so heavy 
top-quark the quark-antiquark annihilation mechanism is dominant over the gluon-gluon fusion 



at Tevatron 



52 



54 



56( 1 and q + b ^ q' + t 



55I ] . So the cross-sections for single top production via q + q' t + b 



53], which are comparable to the quark-antiquark annihilation, are also 



important for the top-quark production. 

On the production of the superhadrons at Tevatron and at LHC, for the reason precisely 
pointed out in the above section we have highlighted the two mechanisms via the hard subpro- 
cesses: gluon-gluon fusion and light quark-antiquark annihilation so far. In fact, there may be 
some other mechanisms for producing the superhadrons which may contribute greater than that 
via light quark-antiquark annihilation and even so sizable to be comparable with that via gluon- 
gluon fusion. For instance, when a comparatively light chargino {'>tl~± < 0{TeV)) is allowed in 
the same SUSY models, the 'single top-squark production' such as that via g + b + X1/2 "^^^ 
be the case: its contribution may be greater than that via light quark-antiquark annihilation 



and even so sizable to be comparable with that via gluon-gluon fusion. It is very similar to 



the top production 



48 



49 



50 



511 ]: at the LHC the single top production via the process 



g + b W + t has a cross section of about 60 pb while the cross section via gluon-gluon fusion 
g + g~^t + t is roughly about 760 pb and the cross section via quark-antiquark annihilation 
q + q^t + t is roughly about 40 pb. Moreover if the bottom-squark bi is also comparatively 
light (m^ < (9(TeV)) in the same SUSY models, then the production via q + q' ti + hi can 
be quite great too. However, all of the possibilities depends on the parameters of the relevant 
SUSY models, thus we would not calculate them precisely in the paper. As for the production 
via the subprocesses such as the annihilation of top-quark and anti-top-quark through gluino g 
or photino 7 exchanging t + i ti+ti, and a top-quark 'scattering' on a gluon: g+t — > ^1 + ^(7) 
etc, we are sure that their contributions to the superhadron production are very tiny due to 
the smallness of the PDF of top-quark in the colliding hadrons. Anyway, for the accuracy 
of the present estimate and such a 'light' top-squark mj^ = 120 ~ ISOGeV, the contribution 
via the quark-antiquark annihilation hard subprocess to the production can be negligible in 
comparison to the dominant gluon-gluon fusion mechanism at Tevatron and LHC both, that is 
quite different from the top-quark case. 

Since the decay of the heavy superhadrons H i.e. (iiQ) with Q = c,b is via the light top- 
squark ti or via the involved heavy quark Q with proper relative decay possibility, so there are 
two typical decay channels for H: one is the decay of the light top-squark ti with the heavy 
quark Q acting as a 'spectator' and the other is the decay of the heavy quark Q with the light 
top-squark ti acting as a 'spectator'. The second decay channel may be quite different from that 
of the decay for a light top-squark ti itself, and then it shall present certain characteristics. 
Therefore, we think that in order to observe and identify (discover) the light top squark ti 
experimentally, one may try to gain some advantages via observing the characteristics of the 
heavy superhadrons H decay. 

If the heavy superhadrons are really observed experimentally, it will be a good news not 
only for the relevant SUSY model (s) but also for the QCD-inspired potential model, because it 
will open a fresh field, i.e. the potential model will need to be extended to treat systems with 
binding of a fermion and a scalar boson. 

As it is known, the fragmentation of a heavy quark 6 or c to a double heavy meson rjf, or 



rjc is quite smaller than that of the heavy quark to a 



possibility about 10"^ ~ lO'^ [22, l23|, M, M, M, 



27 



leavy meson B 01 D i.e. with a relative 



281], thus with the same reason one may 



be quite sure that the fragmentation function of top-squark ti to light superhadrons (tig), 
q = u,d,s are much greater than that of top-squark ii to heavy superhadrons (tiQ) Q = c,b. 



Namely, by conjecture, the fragmentation function of top-squark ti to light superhadrons {tiq) 
can be about (10^ ~ 10^) of the top-squark ti to the heavy superhadron (tiQ) one. With such 
an enhancement this large, the light superhadrons may be produced numerously, and then one 
may collect enough events for experimental observation even at Tevatron. However without 
the additional characteristics due to the decay heavy quark 6 or c of the heavy superhadrons, 
the light superhadrons may be comparatively difficult to be identified experimentally. 

Finally, we should note here that the computation and discussion in the paper are explicitly 
based on the assumption that the light color-triplet top-squark does exist in certain SUSY 
models. As a matter of fact, the results in the present paper are true for a variety of the SUSY 
models, in which even the light top-squark ti is not the lightest SUSY object in nontrivial color. 
As long as in the SUSY models concerned, the hghtest SUSY partner is a scalar in color triplet 
and has a lifetime long enough to form hadrons before decaying, our results as presented here 
remain meaningful by simply replacing the light top-squark ti with the corresponding lightest 
SUSY partner. 
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